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ABSTRACT 

Active galactic nuclei (AGN) with jets seen at small viewing angles are the most luminous and abundant 
objects in the 7-ray sky. AGN with jets misaligned along the line-of-sight appear fainter in the sky, but are 
more numerous than the brighter blazars. We calculate the diffuse 7-ray emission due to the population of 
misaligned AGN (MAGN) unresolved by the Large Area Telescope (LAT) on the Fermi Gamma-ray Space 
Telescope (Fermi). A correlation between the 7-ray luminosity and the radio-core luminosity is established 
and demonstrated to be physical by statistical tests, as well as compatible with upper limits based on Fermi- 
LAT data for a large sample of radio-loud MAGN. We constrain the derived 7-ray luminosity function by 
means of the source count distribution of the radio galaxies (RGs) detected by the ferm/-LAT. We finally 
calculate the diffuse 7-ray flux due to the whole MAGN population. Our results demonstrate that the MAGN 
can contribute from 10% up to nearly the entire measured Isotropic Gamma-Ray Background (IGRB). We 
evaluate a theoretical uncertainty on the flux of almost an order of magnitude. 



1. INTRODUCTION 

The Ferm/-LAT has measured the Isotropic Gamma-Ray 
Background ( IGRB) with very go od accuracy from 200 MeV 
to 100 GeV (lAbdoetal.ll2010cl) . Increased statistics from 
additional LAT data are expected to allow the IGRB to be 
measured over an even broader energy range. The nature of 
the IGRB is still an open problem in astrophysics. Blazars 
and star-forming galaxies contribute 20%-3 0% of the IGRB 
jAckermann et al.l2012aHAbdo et al.l2010dl) . and with a com- 
patible spectral slope. Blazars are Active Galactic Nuclei 
(AGN) whose jets are oriented along the lines-of-sight (l.o.s.). 
Their luminosity is quite high, due to Doppler boosting ef- 
fects. For the same reasons, AGN with axes misaligned with 
respect to the line-of-sight (hereafter MAGN) have weaker 
luminosities but are expected to be more n umerous by a fac- 
tor 2 rr (where Tl is the Lorentz factor) (lUrrv & Padovanil 
119951) . About 10% of the observed AGN are radio-loud. 
According to the unified model, AGN are classified as a 
function of their jet or ientation with respect to the l.o.s. 
JUrrv & Padovanil 1 1995h . A jet misalignment of about 14° 
indicates the separation between blazars and non-blazar, i.e. 
misaligned, objects. In the unified model, radio galaxies 
(RGs) are those objects which, on average, have their jets 
pointing at > 44° from our l.o.s., while MAGN below this 
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angle an d above about 14° are generally identified with radio 
quasars (lBarthelll989l) . RGs are classified in turn into two cat- 
egorie s based upon their radio morphology (iFanaroff & Rilevl 
119741) . The first class of objects, named Fanaroff-Riley Type I 
(FRI), is preferentially found in rich clusters and hosted by 
weak-lined galaxies, and has a low luminosity radio emis- 
sion (which peaks near the center of the AGN and shows two- 
sided jets dominated emission). Fanaroff-Riley Type II (FRII) 
galaxies present a high luminosity radio emission dominated 
by the lobes, while jets and core, when detected, are faint. 
The hot spots, generally not present in FRIs, are usually de- 
tected at the ends of the lobes. The threshold in luminosity for 
discriminating bet ween FRI and FRII is ab out 10^^ W Hz~^ 
sr-i at 178 MHz (IFanaroff & RilevllIgTl . FRIs and FRIIs 
are considered the parent populations of BL Lacertae objects 
(BL Lacs) and flat spectrum radio quasars (FSRQs), respec- 
tively. 

Given the large numbers of known radio-loud MAGN, and 
in some analogy with blazars, it is expected that a non- 
negligible contribution to the IGRB might be attributable to 
the MAGN population. The aim of this paper is the estima- 
tion of the 7-ray diffuse emission produced by the cosmologi- 
cal population of MAGN. We investigate the absolute level of 
the MAGN 7-ray flux and quantify the possible uncertainties 
of our prediction. 

The paper is organized as follows: in Sect. |2]we derive the 
correlation between radio core (Lr,core) and 7-ray luminosi- 
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ties (L^) from a sample of 12 MAGN detected by Fermi-LAI. 
The robustness of the L^~Lr. core correlation is tested in Sect. 
[3] by computing the 95% confidence level (C.L.) upper lim- 
its on the 7-ray flux of a sample of radio-loud MAGN with 
4-years of Fermi data. In Sect. |4]two statistical tests are per- 
formed on the core radio and 7-ray data in order to exclude 
spurious effects in the correlation between luminosities. By 
assuming the relation found between L^ and ir,corc, in Sect. 
|5]we model the 7-ray luminosity function (GLF) from the ra- 
dio luminosity function (RLF). We discuss the consistency of 
the models in Sect. |6l where we compare our predictions of 
the source count distribution and compare them to the Fenni- 
LAT data. Our findings for the contribution of an unresolved 
population of MAGN to the IGRB are presented in Sect. |7] 
together with the evaluation of the relevant uncertainties. Fi- 
nally, we draw our conclusions in Sect. |8] 
Throughout the paper we adopt a standard ACDM cosmology 
with parameters: Hq = 70 km s^^ Mpc^^, ^Im = 0.27, JIa = 
0.73. 

2. THE CORRELATION BETWEEN 7-RAY AND RADIO 
LUMINOSITY 

The calculation of the diffuse emission from unresolved 
{i.e. not detected by the Fer»z/-LAT) MAGN relies on the 
7-ray luminosity function for that specific population. The 
physical processes underlying the emission of 7 rays in RGs 
are not firmly established. It is believed that the bulk of the ra- 
diation is generated via synchrotron self-Compton (SSC) scat- 
terings, where the seed photons are provided by synchrotron 
emiss ion by the same electron population (iMaraschi et aP 
II992I) . It is not excluded that an external inverse Comp- 
ton (EC) scattering occurs o ff photons external to the jet 
jDermer & SchlickeiserllI993l) . Typically the EC scenario re- 
quires that the inverse-Compton (IC) 7 radiation originates 
relatively close to the center of the AGN, within the cen- 
tral parsec, while the SSC 7 rays come from the jet outside 
that region. A relevant contribution to a detectable 7-ray 
emission from the external lobes is t hought to be unlikely . 
The lobe contribution was evaluated in lStawarz et al.l (120061) : 
iMassaro & AJello (2011). Ultra-relativistic electrons in the 
lobes emit synchrotron radiation in the radio band and are able 
to up-scatter low energy photons via IC scattering to high en- 
ergies, provided a high enough electron density is available. 
The dominant contribution is expected to be from CMB pho- 
tons. The IC/CMB scattered emission in the lobes of dis- 
tant galaxies is generally well observed in the X-ray band. 
Extended 7-ray emission from regions correspo nding to ra- 
dio lob es has been detected from Centaurus A (lAbdo et al.l 
l2010al) . Such emission, if interpreted in terms of IC scatter- 
ing of electrons with ambient photons, requires high-energy 
electrons in the lobes, but it is unclear how common this is 
in other RGs. In what follows, we assume the 7-ray radiation 
originates in the central region of the source, as is predicted 
from both SSC and EC scenarios. 

The FRI and FRII galaxies show a strong emission in a wide 
radio band, spanning from hundreds of MHz up to tens of 
GHz. These photons are ascribed to the synchrotron emis- 
sion of highly relativistic electrons moving in the entire re- 
gion of the source. The total radio flux has been measured 
for hundreds of FRI and FRII galaxies. For a number of these 
galaxies the emission from the central unresolved region of 
an arcsecond scale, often referred to as the core, is detected as 
well. In the first and the second catalogs of LAT AGN sources 
dAbdo et al.. ,2010b.: ,Ackermann et al.,,201 lb.) Fermi-LAI has 



reported the detection of 15 MAGN, which can be classified 
into 10 FRI and 5 FRII galaxies (although with some caveats, 
see below). lAbdo et alT ( l2010bl) report on the observation of 
3C 78, PKS 0625-35, 3C 207, 3C 274, Centaurus A, NGC 
6251, 3C 380, 3C 120, 3C 111, 3C 84, PKS 0943-76, while 
Centaurus B, Fornax A and IC 310 have been repo rted in the 
second LAT cata log (2FGL) dNolan et al.ll2012aibl) (for Cen- 
taurus B see also lKatsuta et alj 1 20I3I)) and a Pictor A identi- 
fication has been discussed in lBrown & Adamsl (120 1 2|) . 

In the absence of predictions for the 7-ray luminosity 
function, we follow a phenomenological approach to re- 
late the 7-ray luminosity to the radio luminosity. The lat- 
ter is phenomenologically much better established, given 
the high number of detected MAGN in the radio frequen- 
cies. A possible correlation between radio and 7-ray lu- 
minosities has been proposed for blazars using the En- 
ergetic Gamma Ray Experiment Telesc ope (EGRET) data 
(JPadoyaniet al .11993: Stecker et al. 199 3i:ISalamon & Steckd 
1 1 9941: iDond i & Ghisellini 1995; Naru moto & Totanil 120061) 
Recently, the connection between radio and 7-ray fluxes has 
been explored for both the FSRQs and BL Lacs detected by 
Fermi- LAil during its first year of operation (lAckermann et al.l 
l20Ilai: Ghirlandaet al. 2011). On a similar basis, the rela- 
tion between radio emission and 7-ray data has been studied 
for th ree FRI galaxies observed by EGRET (iGhisellini et al.l 
120051) . as well as for FRI and FRII galaxies with 15 nionths o f 
data taken with Fermi-LAT dlnouelfeoi ItlAbdo et alj|2010bl) . 
Variability studies for FRI galaxies sup port the hypothesis 
of the compactness of the 7-ray source (lAbdo et al.ll20I0bt 
iGrandi et al.l l20I2ah . even if a non-negligible 7-ray coun- 
terpart in radio lo bes has been observed in Centaurus A 
dAbdo et al.ll20I0ah . The situation for the FRII population is 
less de finite. A recent Ferm/ -LAT analysis of the FRII 3C 1 1 1 
galaxy dGrandi et al.ll20I2bl) . together with a multi-frequency 
campaign conducted in the same period, localizes the GeV 
photons from 3C 1 1 1 in a compact, central region associated 
with the radio core. 

The main radio and 7 parameters of all the MAGN observed 
by Fermi-LACr are reported in Table [T] The radio data have 
been chosen to be the closest in time to Fermi-LAI data tak- 
ing. Whenever a significant variability has been found, we 
have selected radio data almost contemporary to the 7-ray ob- 
servations. Radio data have been taken with the Very Large 
Array (VLA) for all the objects except NGC 6251, measured 
with the Very Long Baseline Interferometer (VLBI). The lin- 
ear size scales explored by the instruments depend on the red- 
shift of the sources. In our sample, it varies from about 0.01 
kpc to a few kpc, except for NGC 6251 (0.002 kpc) and 3C 
380. Data for 3C 380 are taken from Effelsberg observations 
with a high resolution. However, this source shows a compact 
steep spectrum radio morphology and the radio flux from the 
central region is close to the total emission and to the flux 
measured with a few arc sec scale Q. For 3C 84 the variabil- 
ity is very pronounced and we have therefore excluded it from 
our correlation analysis. IC 310 lacks measurements of the 
core at 5 GHz, and the total radio flux is very faint. For PKS 
0943-76 only upper limits for the core are given. For these 
reasons, these two galaxies are listed but not considered in our 
analysis. The photon index F valid between 0.1 and 10 GeV, 
and the fer»z/-LAT flux integrated for E^ >0.1 GeV have 
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been taken fr omlAbdoet al.l (i20I0bl) for 3C 78, 3C 111, 3C 
120 and from lNolan et al.l d2012bh for the remaining objects. 



http://3crr.extragalactic.info/cgi/database 
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Table 1 

Main radio and 7 properties of the MAGN observed by Fermi-hAI. Column 1 : name of the MAGN (radio classification: FRl or FRIl). 2: redshift, 3: Galactic latitude, 4: spectral index for radio core (total) 
spectrum in a range including 5 GHz, 5: measured radio core (total) flux at 5 GHz; 6: photon index for 7-ray spectrum between 100 MeV - 100 GeV; 7-7-ray flux above 100 MeV; 8- Radio core luminosity at 5 

GHz ; 9- 7-ray luminosity 

Referen ces: iM organti et al. ( 1993): 2 Pauhnv-Toth et al. 1 1972); 3 Ekers et al. ( 1989); 4 Kuehr et al. 1 1981); 5 Spinrad et al. 1 1985); 6 Muflin et al. (2006); 7 Naear et al. (2001); 8 Giovannini et all (7981) : 
9 <Laing et al. ( 1983); 10 Israel et al. (2008); UBums et al. (1983); 12-Wright et al. ( 1994); 13 Evans et al. (2005); 14Mantovani et al. (2009); 15 Jones et al. (2001); 16-Massardi et al. (200l; 
17 4Geldzahler & FomalonJ 119841) ; IS jBecker et al.l ( ll991l) ; 19 jPerlev et al.H1997.') : 20 Linfield & Perley ( 1984); 2 UKadler et alj ( 120121) ; 22 jBurgess & HunsteadI ( 120061) ; 23 jKadler et alj )20'Ta ; 

24 .Gregorv & Condon C199U 



MAGN(FRI,FR1I) 


z 


bn 


acoro(Qtot) 


S^^i^?^ [Jy] (S^^i^n [Jy]) 


r 


F^ [10"9phcm-2s-i] 


L5c>H^,[ergs-i] 


L^[ergs-i] 


3C78/NGC 1218(1) 


0.0287 


-44.6 


0(0.64 1) 


0.964 ±0.1641 (3.40 ± 0.112) 


1.95 ±0.14 


4.7 ±1.8 


(8.8 ± 1.4 ) ■ lO'"' 


(1.11 ±0.54) ■ 10-13 


3C 274/M 87(1) 


0.0038 


74.5 


0(0.79**) 


3.0971 ± O.OSOO'^ (71.566 ± 0.993 ^) 


2.17 ±0.07 


25.8 ±3.5 


(4.90 ± 0.05 ) ■ 10^9 


(6.2± 1.1 ) ■ 10*1 


Cen A(l) 


0.0009 


19.4 


0.3010 (0.70^0) 


6.984 ± 0.210" (62.837 ± 0.099 i^) 


2.76 ±0.05 


175 ± 10 


(6.19 ±0.19) ■ lO^* 


(1.14 ±0.09) -10*1 


NGC 6251(1) 


0.0247 


31.2 


0(0.729) 


0.38 ± 0.04 13 (0.510 ± 0.050 i^) 


2.20 ±0.07 


18.2 ±2.6 


(2.57 ± 0.27 ) ■ 10-10 


(1.82 ±0.41 )-10*3 


Cen B(l) 


0.0129 


1.68 


0(0.13i«) 


2.73015 (6.58 ±1.0416) 


2.33 ±0.12 


39.3 ±11.4 


5.02- 10-10 


(8.6 ± 3.2 ) ■ 10*2 


For A(I) 


0.00587 


-56.7 


0.501^(0.52 1) 


0.051" (72i) 


2.16 ±0.15 


7.7 ±2.4 


1.93 10^* 


(4.6 ± 2.2) • 10*1 


3C 120(1) 


0.0330 


-27.4 


0(0.441*) 


3.458 ±0.588 1 (8.60±1.46i) 


2.71 ±0.35 


29 ±17 


(4.20 ±0.71 ) ■ 10*1 


(2.9 ± 1.6) ■ 10*3 


PKS 0625-35(1^ 


0.0546 


-20.0 


0(0.65-') 


0.600 ±0.030 3 (2.25 ± 0.09 '') 


1.93 ±0.09 


12.9 ±2.6 


(2.02 ±0.10) ■ 10-11 


(1.21 ±0.43) ■ 10** 


Pictor A(II) 


0.0351 


-34.6 


0(1.07 1) 


1.15 ±0.0519 (15.45 ±0.47") 


2.93 ±0.03 


21.9 ±3.6 


(1.58 ±0.07) ■ 10-11 


(2.13 ±0.46) ■ 10*3 


3C 111(11) 


0.0485 


-8.61 


-0.20^(0.73^) 


1.1420 (6.637 ±0.9961*) 


2.54 ±0.19 


40 ±8 


2.98-10*1 


(1.01 ±0.38) -10** 


3C 207(11) J^ 


0.681 


30.1 


0(0.90"^) 


0.5391 ± 0.0030*' (1.35 ± 0.04 ") 


2.36 ±0.11 


17.3 ±3.3 


(3.32 ± 0.02 ) ■ 10-1=* 


(2.41 ± 0.61 ) ■ 10**5 


3C 380(ll)t 


0.692 


23.5 


0(0.71^) 


5.073 ± 0.105 " (7.45 ± 0.37 *) 


2.34 ±0.07 


30.3 ±3.7 


(3.12 ±0.07) ■ 10-1-1 


(4.44 ± 0.73 ) ■ 10**5 


IC 310(1) 


0.0189 


-13.7 


n.a.(0.752^) 


n.a. (0.258 ±0.0312'!) 


2.10 ±0.19 


11.1 ±6.2 


- 


(7.9 ± 4.9 ) ■ 10*2 


3C 84/NGC 1275(1) 


0.0176 


-13.2 


(0.785) 


high variability 


2.00 ±0.02 


175 ±8 


- 


(1.22 ± 0.07 ) ■ 10** 


PKS 0943-76(11) 


0.270 


-17.2 


n.a. 


upper limits(0.75722) 


2.44 ±0.14 


19.5 ±5.1 


- 


(2.47 ±0.71 )-10*5 
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Figure 2. Upper limits on several Ferair'-LAT undetected radio-loud MAGN. 
Blue (red) symbols refer to possible FRI (FRII) classification. The correlation 
in Eq. |5](solid black line) is displayed together with the la eiTor band (light 
pink shaded area). 

have calculated luminosities according to Eqs. [T]|4] propagat- 
ing errors on F while neglecting errors on the redshift, given 
their negligible effect . The correlation between ir,corc and 
L^ for the 12 objects is described by the function: 



SGHzN 
r.corc/ 



(5) 



log(L^) = 2.00±0.98 + (1.008±0.025)log(L 

represented by the solid line in Fig. [T] while the relevant la 
error band is shown as a shaded area. It is obtained from the 
errors on both the 7-ray and radio luminosities, but the former 
dominates the uncertainty. Eq. |5]describes a linear correlation 
in the log-plane with a coefficient very close to one. 7-ray 
luminosities are greater than radio luminosities by about two 
orders of magnitude. The uncertainty band of the 7-ray fluxes 
measured by the Fermi-LNI spans one order of magnitude 
around the best fit. The significance of the correlation is tested 
in SecH 

It is useful to compare the Lr,coro - L-y correlation obtained 
by removing three sources with non standard properties from 
the set of 12 MAGNs in Tabl^D The radio morphology of 3C 
380 is not unambiguously determined, as it shows some prop- 
erties of a compact steep spe ctrum radio source (lAbdo et al.l 
I2010bt lWiIldnsonetaI]ll99l . The FRII 3C 207 behaves as 
a steep spectrum radio quasar in the optical band (lAbdo et al.l 
I2010bl) . while PKS 0625-35 has no clear association. We 
have therefore also calculated the correlation function exclud- 
ing the galaxies 3C 380, 3C 207 and PKS 0625-35: 



log {L^) = 2.1 ± 2.1 + (1.005 ± 0.055) log {LI 



5GHzN 



(6) 



The result is not very different from Eq. |5] if it were not for 
the increased spread in the fitted coefficients. Indeed, 3C 207 
and PKS 0625—35 have large errors and 3C 380 is quite close 
to the correlation with the whole sample. 

We report here also the correlation between the total radio 
luminosity at 5 GHz and the 7-ray luminosity for the whole 
sample of 12 sources. The correlation is found to be: 



log {L^ 



-2.5±1.1 



(1.095 ±0.026) log (i5GHz)_ (7) 



The experimental values for the total radio luminosity are 
quoted in Table [T] The fit for the sample of 9 sources results 
with: 



9 sources T he correlation impUed by Eq.|5]is close to the one 
obtained in iGhisellini et al.l (l2005h for the very small sample 
of EGRET 7-ray loud FRI galaxies. In the case of blazars 
the slope of the correlation between L-y(> 100 MeV) and 
radio luminosity at diffe rent frequencies was found to be: 
1.07± 0.0 5 at 20 GHz dGhirlanda et al.1 I20T0I). 1.2 ± 0.1 
at 5 GHz (ISteckeretal.l ll993b and 1.06 ± 0.02 at 8.4 GHz 
dYuan & Wang|l20I2b . The slope coefficient of the correlation 
for RGs is therefore similar to the correlation for blazars. This 
might indicate that the 7-ray emission mechanism is similar 
for MAGN and blazars. We therefore assume that the correla- 
tion in Eq. |5]is a good representation of the luminosity of the 
cores of MAGN and we will employ it in the remainder of this 
work in order to derive the emission of the MAGN population 
not detected by the Fenni-l^P^, but potentially providing a 
non-negligible diffuse flux. 

3. UPPER LIMITS FROM RADIO-LOUD FRI AND FRII GALAXIES 
NOT DETECTED BY FERMI-hAT 

In order to test the robustness of the core radio- 7 correla- 
tion found in Eq. |5]we study a sample of radio-loud FRI and 
FRII galaxies that have not been detected by Ferm/-LAT. For 
these objects we derive 95% C.L. 7-ray upper limits and ver- 
ify that they are consistent with Eq. |5] given the uncertainty 
band s hown in Fig. [T1 The sar nple has been extracted from 
RGs in lKataoka et alI ( l20Tlh and lGhisellini et alj ( l2005t) (first 
and second block in Table |2l respectively), and represents the 
sources with the highest radio core fluxes at 5 GHz. Further 
selection criteria have been applied in defining the sample for 
our purposes. From the sample of broad line RGs whose up- 
per limits have been pr esented in iKataoka et al.l (1201 Ih we 
have excluded Pictor A (iBrown & Adamsll2012l) . detected in 
the meanwhile, and the sources that do not show a clear FRI 
or FRII radio morphology classification (RGB J1722H-246 
and PKS 2251H-11 being Seyfert galaxies, S5 2116H-81 be- 
ing a flat spectrum radio source with a radio jet morphol- 
ogy). Moreover, sources with latitudes below 10 ° have been 
rejected in order to avoid a strong contamination from the 
Galactic plane foreground. This criterion applies to 4C 50.55 
(b = 0.39°). Th e same criteria have been applied to sources in 
iGhiseUini eFaP (12005 ) leading to the exclusion of 3C 84, 3C 
274, 3C 78 already detected in 7 rays, and 3C 75 that has an 
atypical RG morphology . Finally, 3C 317 ha s been excluded 
because of its variability (IVenturi et al.l2004l) . Four FRII RGs 
from the 3CRR catalogue (3C 245, 3C 109, 3C 212, DA 240) 
have been added to the sample in order to cover a wider range 
in radio luminosity (last block in Table |2]i. Our sample is 
therefore composed of 17 FRII and 22 FRI RGs. 
We have computed 7-ray flux upper limits for the listed galax- 
ies by using the Fenni-LAT Science Tools Q. The data taking 
period for the analysis is from the starting time of the mission, 
2008 August 4, until 2012 September 9. The Mission Elapsed 
Time (MET) interval runs from 239557414 to 368928003. 
Data have been extracted from a region of interest (ROI) of 
radius = 8 ° centered at the position of the source. This radius 
represents the best angular region for sou rce analysis as long 
as sources are far from the Galactic plane (lAbdo et al.l2009bF . 
and indeed we neglect in this analysis sources that lie below 

10 ° in latitude. We selected 7-rays in the energy range 100 
MeV - 100 GeV. 



log {L^ 



3.5 ±2.3 



(0.948 ±0.056) log (L^™^)- (8) 



' http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation, 
version v9r27pl, Instrumental Response Functions (IRFs) P7_V6 



software 
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We are using P7SOURCE_V6 photons. Good survey data are 
selected accordingly to software recommendations, with the 
rocking angle selected to be less than 52°. Data selection and 
preparation eliminate photons from the Earth limb by apply- 
ing a cut on the zenith angle of 100 °. An unbinned maximum- 
likelihood analysis was performed. In the cases where the fit 
did not converge we have performed a binned analysis as rec- 
ommended. We therefore analyze the source region with both 
methods and draw the upper limits with the help of the LA- 
TAnalysisScripts 0, which make use of the UpperLimits.py 
module. 

Each galaxy in the sample was modeled as a point-like source 
with a power-law spectrum of index F = 2.5. This value has 
been cho sen as nominal spectr al index for all MAGN in anal- 
ogy with iKataoka et all ( 120111) . We have verified that choos- 
ing r = 2.3, closer to the distribution of the spectral indices 
from [T] changes the limits by ^ 10%, while an index of 2.7 
leaves results unchanged. The number of expected counts 
in the ROI is derived by considering the emission from all 
sources in the 2FGL inside a source region (distance from 
the target region) of 13° (8° + 5°). The fitting procedure 
leaves the spectral parameters of all the sources inside the 
ROI free, whereas sources in the region 8° < r < 13° have 
spectral parameters fixed to the values of the 2FGL. Addi- 
tional backgrounds are the Galactic diffuse emission and the 
isotropic diffuse model, which includes the true IGRB and the 
residual particle contamination Q The diffuse models used in 
the analysis are: gal_2yearp7v6_v0.fits for the Galactic dif- 
fuse model and iso_p7v6source.txt for the isotropic spectral 
template. All the relevant normalizations have been left as 
free parameters during the fitting procedure. The method used 
to compute the upper limits is a standard profile likelihood. 
A 95% upper limit (UL) has been computed when the Test 
Statistic (TS) was less than 25. In Tabl^the flux upper lim- 
its are quoted together with the TS value for both unbinned 
and binned analysis. Given the systematic uncertainty arising 
from the different statistical methods, we adopt as upper limit 
the highest value for the flux bound. These conservative upper 
limits are shown in Fig|2]with the luminosity correlation from 
Eq. |5] overlaid. It is evident that the calculated upper limits 
do not fall below the uncertainty band, thus corroborating our 
core radio and 7-ray correlation. 

4. TEST OF THE RADIO-7 CORRELATION 

The correlation established in Eq.|5]could be biased by dis- 
tance dependence of the luminosity and flux-limited samples. 
We have tested the strength of the correlation via a partial cor- 
relation analysis, in order to verify that the radio core - 7-ray 
luminosity correlation for MAGN is not spurious. 

Firstly, we calculate the Spearman rank-order correla- 
tion coefficient. The Spearman correlation coefficients are 
0.94, 0.92, 0.98 between log (i^™|) and log {L^), between 
log (i^™^) and redshift, and between log (L^) and redshift, 
respectively. The partial correlation coefficient turns out to 
be 0.51 and the null hypothesis that the two luminosities are 
uncorrected is rejected at the 95% C.L. 

Moreover, we test the significance of the radio-7 correlation 
by using the modified Kenda ll r rank correlation test proposed 
by lAkritas & Siebeij (Il996h . which is suitable for partially - 

^ User contributions http://fermi.gsfc.nasa.gov/ssc/data/analysis/user/ 

' http://fermi.gsfc.nasa.gov/ssc/data/access/lat/2yr_catalog/ 

* http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 
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Kendall T Correlation Coelficient 

Figure 3. Null hypothesis distribution of t correlation coefficients assum- 
ing independence between the 7 and radio wavebands. The null hypothe- 
sis distributions are generated from 5800 pemiutations of 7-ray luminosities 
among the galaxies by requiring that the resultant 7-ray fluxes exceed the 
flux threshold of 7.7 10^^ photons cm^^ s^^. The mean value is 0.223 
with the standard deviation of the distribution RMS=0. 173. The correlation 
coefficient of the actual data is represented by the red solid line, r = 0.397. 

ce nsored datasets. By foll owing the procedure highlighted 
in lAckermann et aH (l2012al) . we perform a Monte Carlo sim- 
ulation in order to compute the distribution of the r corre- 
lation coefficients obtained under the null hypothesis of in- 
dependence between the two wavebands. Starting from the 
detected sample of 12 sources we build several dataset real- 
izations by scrambling the derived 7-ray luminosities among 
galaxies. For each galaxy we then compute the correspond- 
ing flux and we retain only galaxies with a flux above the 
minimal 7-ray flux of the detected sample (7.7 -10"^ photons 
cm^-^ s^^). If the scrambled sample has fewer than 12 sources 
above the flux threshold, we randomly extract an additional 
source from the upper limit dataset (from Table |2) until the 
flux threshold is reached. For each scrambled dataset we then 
compute the Kendall coefficient and we build its distribution 
as shown in Fig. [5] The displayed distribution refers to 5800 
realizations of scrambled samples and the red line represents 
the value of the r correlation coefficient of the actual data, r 
= 0.397. 

Finally, we compare the t correlation coefficient of the ac- 
tual data to the distribution of r and we find that the integral 
of the distribution above r = 0.397 is 0.05. This is the prob- 
ability to obtain the actual correlation by chance, i. e. the 
p-value of the correlation (the smaller the p-value, the greater 
the probability for the observed correlation of being true). As 
in the case of the Spearman test, we can exclude the correla- 
tion happening by chance at the 95% C. L . The result indi- 
cates a physical correlation between the core radio emission 
and the 7-ray flux of the Ferm/-LAT detected MAGN. 

5. THE 7-RAY LUMINOSITY FUNCTION 

The luminosity function for a given energy defines the num- 
ber of sources emitting at that energy per unit comoving vol- 
ume, per unit (base 10) logarithm of luminosity: 



piL,: 



d^N 



dlog{L) dV' 



(9) 



In the radio band, data are available for hundreds of radio- 
loud MAGN, depending on the frequency of the survey. Usu- 
ally radio observations refer to the total emission of the AGN, 
including the central region, jets and radio lobes. Only for a 
limited number of objects detected at low radio frequencies 
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Table 2 

Flux upper limits on a sample of MAGN. Column 1 : name of the MAGN (radio classification: FRI or FRII), 2: redshift, 3: measured radio core flux at 5 GHz 

[Jy], 4: TS of unbinned analysis, 5: 95% C.L. upper limit from unbinned analysis on the flux above 100 MeV in units of 10^^ ph cm^'^ s^^, 6: TS of binned 

analysis, 7: 95% C.L. upper limit from binned analysis on the flux above 100 MeV in units of 10^® ph cm~^ s~^; 8: radio core luminosity at 5 GHz in units of 

ergs~^. 

References: UMorganti eFanfT^^ : 2 Wenstock eTZIH^g^: 3-Third Camb rid ge Catalogue of Rad io Sources ; 4 jDodson et aTH^OnSI) : 5 4Neff et aljfT995l) : 

6 4Pearson etal] n992h : T JTingav et all (imS ). 



MAGN(FRI,FRII) 


z 


^icHz [Jy] 


'-'^unbinned 


pUL 

unbinned 


TSbi7ined 


pUL 

binned 


rSOHz 
^r.covc 1 


3C 18 (II) 


0.188 


0.083^ 


< 1 


2.7 


2.6 


6.0 


3.51 


10-ii 


B3 0309+41 IB (II) 


0.134 


0.320^ 


- 


- 


< 1 


5.8 


6.73 


10-ii 


3C215(II) 


0.412 


0.0164-^ 


<i 


3.1 


4.1 


6.0 


3.56 


lO^i 


3C 227 (II) 


0.086 


0.032^ 


< 1 


0.1 


< 1 


1.1 


2.70 


IQJO 


3C 303 (II) 


0.141 


0.1 50^^ 


< 1 


2.8 


3.3 


4.6 


3.50 


lO^i 


3C 382 (II) 


0.058 


0.1 88^* 


< 1 


4.1 


1.2 


5.9 


7.12 


1040 


3C 390.3 (II) 


0.056 


0.1 20^ 


< 1 


1.7 


3.0 


4.7 


4.26 


10^" 


3C 411(11) 


0.467 


0.078'5 


- 


- 


< 1 


6.1 


2.2- 


10^^ 


4C 74.26 (II) 


0.104 


0.1 00^ 


1.1 


5.4 


<1 


5.7 


1.25 


lO^i 


PKS 2153-69 (II) 


0.028 


0.300 ^ 


4.2 


6.6 


<1 


6.2 


2.67 


1040 


3C 445 (II) 


0.056 


0.086^ 


<1 


0.8 


<1 


1.0 


3.06 


1040 


3C 465 (I) 


0.029 


0.270» 


- 


- 


< 1 


0.5 


2.5- 


lO^u 


3C 346 (I) 


0.162 


0.220^ 


4.5 


6.4 


10.8 


10.2 


1.39 


1039 


3C 264 (I) 


0.021 


0.200^ 


9.0 


5.7 


14.0 


7.5 


9.58 


1039 


3C 66B (I) 


0.022 


0.1 82^^ 


- 


- 


< 1 


8.3 


9.31 


1039 


3C 272.1(1) 


0.003 


0.180-^ 


5.2 


5.6 


5.3 


6.8 


1.66 


1038 


3C315(I) 


0.1083 


0.1 50^^ 


- 


- 


< 1 


2.1 


2.04 


10-li 


3C 338 (I) 


0.030 


0.105-^ 


- 


- 


< 1 


4.6 


1.07 


104u 


3C 293 (I) 


0.045 


0.100^ 


< 1 


1.5 


<1 


1.8 


2.29 


104u 


3C 29 (I) 


0.045 


0.093-' 


< 1 


1.5 


<1 


4.1 


2.11 


1040 


3C31(I) 


0.017 


0.092^ 


- 


- 


< 1 


4.0 


2.83 


1039 


3C310(I) 


0.054 


0.080-' 


< 1 


1.2 


<1 


2.1 


2.63 


104U 


3C 296 (I) 


0.024 


0.077-' 


< 1 


1.5 


<1 


2.3 


4.79 


1039 


3C 89 (I) 


0.1386 


0.049-' 


- 


- 


<1 


1.8 


1.10 


lO^i 


3C 449 (I) 


0.017 


0.037-' 


< 1 


0.5 


<1 


0.8 


1.19 


1039 


3C 288 (I) 


0.246 


0.030-' 


< 1 


1.5 


1.6 


3.7 


2.22 


lO^i 


3C 305 (I) 


0.0414 


0.0295'^ 


- 


- 


<1 


2.1 


5.66 


1039 


3C83.1B(I) 


0.026 


0.040-* 


10.0 


19.7 


16.5 


23.2 


2.89 


1039 


3C 424 (I) 


0.1270 


0.0180 ^ 


- 


- 


<1 


1.6 


3.39 


10* 


3C 438 (II) 


0.290 


0.007 1-' 


< 1 


0.9 


<1 


3.2 


7.40 


1040 


3C 386 (I) 


0.018 


0.120-^ 


- 


- 


<1 


3.2 


4.15 


1039 


3C 277.3 (I) 


0.0857 


0.0122-^ 


- 


- 


4.2 


5.1 


1.03 


104U 


3C 348 (I) 


0.1540 


O.OIO^ 


- 


- 


<1 


5.1 


2.80 


1040 


3C 433 (II) 


0.102 


0.005^ 


- 


- 


<I 


1.9 


5.96 


1039 


3C 442A (I) 


0.027 


0.002» 


<1 


0.7 


<1 


0.9 


1.62 


1038 


3C 245 (II) 


1.029 


0.910^* 


< 1 


2.0 


<1 


4.0 


1.30 


1044 


3C 109 (II) 


0.306 


0.263» 


< 1 


1.4 


<1 


3.5 


3.06 


10^^ 


3C212(II) 


1.049 


0.150» 


6.4 


7.1 


10.11 


8.8 


2.22 


10^^ 


da 240 (II) 


0.036 


0.105-^ 


< 1 


1.5 


<1 


2.8 


1.48 


104u 



(around 0.1 -few GHz), the flux from the central core alone 
has been measured. The RLF is derived phenomenologically 
by fitting data on the emission of the radio sample. Results on 
the total RLF are quite well established dWillott et al.l 120011: 
iDunlop & Peacockll 19901:1 Yuan & Wang|l2012h while the lit- 
erature about the core radio luminosity function is still limited 
(I Yuan & Wand20r2b . given the scarcity of experimental data. 
Unfortunately, deriving the GLF from fitting the gamma- 
ray measurements is not feasible, due to the small size 
of the 7 -ray loud MAGN sample. Following previ- 
ous att empts applied to blazars ([Stecker & Salamoi 
Kazanas & Perl mari 11997^ ' 



7-ray luminosities found in Sect. |2] We assume that; 



N^^kNr 



(10) 



1996; 



2006; 



" I2009F 
and, to 



Inoue & Totani 

Stecker & Venters 2 0111) 

^Ghisellini et al...2005nin"ouell201 lb 
the RLF by exploiting the 



ami 
IJ 



where the normalization k takes into account our ignorance 
of the number of radio-loud MAGN emitting in 7 rays as 
well {Nr and N^, respectively). From Eq.|9] it follows that 
N ~ J dV J p{L, z)d log L and therefore the GLF is defined 
through a RLF by: 



Narumoto & Totanil 

Abaz aiian et al. 20 IJ 

a lesser extent, to RGs 

, we derive the GLF from 

correlation between radio and 



p^{L^,z) = k Pr{Lr,z) 



d log Lr 



(11) 



d log L^ 
Given the results of the previous sections, the above equation 
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takes the specific form: 



P^{L^, z) — k Pr,corG(i'r,corc(^7)' ^) 



dlogL5™^,(L^) 



d log L^ 



(12) 

where Pr.coic refers to the radio luminosity function of the 
cores of the MAGN. If our hypothesis of a correlation be- 
tween the core radio and 7 emission is physical, as supported 
by the results on the ULs (see previous section), we might ex- 
pect k values not too far from 1. In other words, each RG 
with a bright radio core is expected to emit in the 7-ray band 
as well. The correlation between radio and 7-ray luminosi- 
ties is assumed to be a specific analytical expression, Eq. |5] 
shown to be in very good agreement with the data. In this 
sense, the scatter in the correlation derives only from errors in 
the experimental data and not in a potential scatter on the lu- 
minosity form. The radio luminosity is energetically weaker, 
according to Eq. |5] As already noted, the lack of a reliable 
core RLE from data prevents us from using EqJTSjdirectlv. We 
will therefore make use of the total RLE and obtain the core 
RLE through the link between total and core radio luminosi- 
ties. 

As a first ingredient, we need a correlation between radio core 
and total luminosities. In Eig|4]we display the correlation be- 
The three curves correspond to: 



tween L^™^ and L?gh-„ 



1.4GHz 



(13) 



r408MHz 



(14) 



(15) 



logL^^^,: = 4.2 ± 2.1 + (0.77 ± 0.08) logL^ 

from lLara et al.l (120041) (black solid curve), 

logL^™/, = 7.6 ± LI + (0.62 ± 0.04) logL^ 

from lGiovannini et al.l (1200 lb (pink dot-dashed curve), 
logL^osMHz ^ 7.io±0.90 + (0.83±0.04)logL5GHz 

from lYuan & WangI (120121) (green dotted curve). We report all 
the luminosities at 5 GHz, assuming a power-law dependence 
L/i/ ex z^^", with atot = 0.80. It is clear from Eig. |4]that 
the experimental data for our MAGN sample are best r epre- 
sented by the correlation proposed bv lLara et al.l (120041) . We 
will therefore adopt Eq{T3]throughout the paper The possible 
systematics introduced by this correlation are likely compen- 
sated, at least to a good extent, by the fit to the cumulative 
number counts (see the following section). 

The core RLE may be obtained from the total one following 
the same reasoning as for Eqs. [TOlfTTI 



Pr,i 



^(-tJr.coroi z) — Pr.tot (-t^r,tot 5 Z) 



d log ir, 
dlogLr^t 



(16) 



where we expect that the number of MAGN showing core and 
total emission is almost the same. We adopt as the total RLE 
the one derived in lWillottet all (120011) (Model C with Oa/=0) 
and shift luminosities from 151 MHz to 5 GHz according to 
the power law explained above. We convert the comoving 
volume to the standard ACDM cosmology by using the con- 
version factor ri{z): 



r]{z) 



(PVw/dzdn 
dW/dzdfl ' 



(17) 



where d^Vw /dzdi l is the comoving volume element used by 
lWillottetalJ(l200l : 



vw 



c'^z^{2 + zy 



dzdn 4ff3^(l + z)3' 



(18) 



5 26 



Fermi FRI ^ 
Fermi FRil r 
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Giovannini 
Yuan and Wang 



Figure 4. Radio core luminosity v ersu s total radio luminosity at 5 GHz. 
Solid black line corresponds t o E g. |13| whi le the pink dot-dashed and the 
green dotted correspond to Eg |14l and Eq llSl respectively. Blue squai'es (red 
open circles) are the experimental data for our sample of FRI (FRII) taken 
from Table [T] 

c is the speed of light and ffo.w = 50 km s~^ Mpc~^. In the 
cosmological model ACDM the comoving volume element is 
defined as: 



d^V 
dzdfl 



(19) 



cdL^z) 



Ho{l + z) V(l - ^A - ^A/)(l + 2)2 + (1 + zfilM + ^A 

We finally obtain the GLE inserting EqlT6]in EqfT2l 

Py{L^,z)=k Pr^tot (ir,to/(-^r,coro(-^7))'^) 
dlogLSGHz dlogL^GHz 



dlogL^ dlogi^GHz 



(20) 



The d log i^™^/(i log L^ will be computed from Eq|5] while 
the (ilogL^™^/dlogL^™.^ derives from the total-core cor- 
relation, EqfT3] 

6. PREDICTIONS FOR THE SOURCE COUNT DISTRIBUTION 

An important observable for the correctness of our method 
is provided by the source count distribution of MAGN mea- 
sured by Fenni-LAil. The source count distribution, known 
also as log A^ — log S, is the cumulative number of sources 
iV(> F^) detected above a threshold flux F^. We have de- 
rived the experimental source c ount distribution o f the 12 
MAGN of our sample following (lAbdo et al.ll2010el) : 



N{>F-,,0 



N{>F,)^ J2 



1 



^{F^.i)'' 



(21) 



where the sum runs on all the i-sources with a 7-ray flux 
F^^i > Fj, and u!{Fy^.i) is the flux dependent d etection ef- 
ficiency compa tible with our sample. As shown in lAbdo et al.l 
(I2009ail2010eh . at faint fluxes the Fenni-LAT more easily de- 
tects hard-spectrum sources rather than sources with a soft 
spectrum. Sources with a photon index of r=1.5 can be de- 
tected down to fluxes that are a factor > 20 fainter than those 
of a source with a photon index of 3.0. Given this strong 
selection effect, the intrinsic photon-index distribution is nec- 
essarily different from the observed one. This effect is taken 
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into account by the detection efficiency. Since the latter is not 
available for the MAGN sample, we re asonably assume it is 
the same as for blazars and take it from lAbdo et al.l(l2010db . 
The theoretical source count distribution A^th(> ^7) for a 7- 
ray flux K^ is calculated following the definition of GLF in 
Eq.|9l 



A^th(>i^.)=47r 



dN 

It 



dV 



dL^ 



L.,{F^,z.T) i7ln(10) 



-- d^V 
dzd^ 

p^{L^,z,T), (22) 



where L^ {F^ , z, F) is the 7-ray luminosity of a RG at redshift 
z, whose photon spectral index is F and photon flux is F-y 
(integrated above 100 MeV). The spectral index distribution, 
dN/dT, is a ssumed to be gaussian in analogy with blazars 
(lAbdoetal.ll2010dl) . The comoving volume, d^V/{dzdVi), is 
computed according to Eq{T9] We fix r,„i„=1.0, F,„aa;=3.5, 
Zmax = 6 and L^^rnax = 10^" erg s~^ 

Fig. [5] shows the theoretical iVth(> ^7), calculated from 
Eql22l with several bands of uncertainty, overlaid with the ex- 
perimental source count distribution from EqJ2T] Their com- 
parison is discussed here only as a consistency check of the 
validity of the assumptions involved in Eq|22]and in particu- 
lar of the ratio of MAGN emitting in 7 rays relative to those 
emitting in radio-core, i.e. the k parameter in Eq. |22] The data 
points for the experimental A^(> F^) are in fact highly corre- 
lated, and a fit to those points is not statistically meaningful. 
Nonetheless, it is useful to fit the theoretical A^tii(> F^) to the 
experimental source count distribution to constrain the only 
free parameter k in this calculation. Additionally, the shape 
of the function predicting A^(> F^) is essentially driven by 
the radio luminosity density function, and not by the fit to the 
experimental source count distribution. 
The black dashed line in Fig. |5]has been derived from the 
best fit parameters of the ir.corc — L^ correlation (namely 
Eq|5]l, whose fit to the experimental source count distribution 
gives k — 3.05. This indicates that the best fit radio core-7 
correlation function slightly under-predicts the distribution of 
MAGN observed by Fermi-LAT. 

For obtaining the bands depicted in Fig. |5] we have proceeded 
as follows: 

i) we have calculated the A^(> F^) for all the correla- 
tion coefficients falling in the la uncertainty band for the 
Lr,corc ~ Lj relationship (Fig. [U; 

m) for each combination of these coefficients we have deter- 
mined k from the comparison with the log N — log S (pink 
shaded area); 

in) the configuration with the lowest x^ among all the config- 
urations explored at point ii predicts the best iV(> F^) (black 
soHd fine, fc=0.258); 

iv) all the configurations giving a Icr variation from the low- 
est x^ (minimal x^ H-3.53) span the cyan shaded area. 
The red dot-dashed curve was obtained for the radio total 
- 7-ray luminosity correlation in Eqi?] and the total RLE in 
IWillott et al.l (120011) . This hypothesis leads to a lower number 
of sources at the lowest fluxes. 

Finally, the green shaded b and wa s obtained by fixing the nor- 
malization factor k in Eqs. llOllll equal to 1, which represents 
the ideal situation in which we predict that each MAGN has 
a radio-loud central region emitting in 7 rays as well. This 
result is an important test of the validity of our initial assump- 
tion that a MAGN with a radio core emission also emits pho- 



tons in the 7-ray energy band, via likely SSC and EC pro- 
cesses. It is remarkable that the band is a good fit to Fermi- 
LAT data. 

Given the uncertain classification of some of the sources, 
as explained in Sec|2] we also provide the source count distri- 
bution for the 9 sources with firm FRI or FRII classification, 
Fig|6] We show the experimental and the theoretical source 
count distribution predicted when the three galaxies 3C 380, 
3C 207 and PKS 0625—35 are excluded from the analysis. 
The black solid line is the same as in Fig. |5] but obtained with 



— L-y lU 



9 data points and employing Eq. |6] for the Lr 
minosity correlation function. The result is compatible with 
data, with x^ =4.65 and the normalization for the source num- 
ber distribution fc=2.37. The red dot-dashed curve is the same 
as in Fig. |5] but obtained from the total RLE and EqlS] and 
minimized with respect to the 9 data points. 

7. THE DIFFUSE 7-RAY EMISSION FROM MAGN 

The diffuse 7-ray flux due to the whole population of 
MAGN may be estimated as follows: 



d^Fje) 
dedfl 



^-^^d_N 
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dzdVL 



dF^ 
IT 



(23) 



L..,ln(10) 



P-yiLj, z)(l - uj{F^{L^,z))) exp {-T^^^{e, z)). 



The minimum 7-ray luminosity value is set to 10^^ erg s^^, 
the maximum at 10^*' erg s^^. The term uj{F^{L^,z)) is 
the detection efficiency of Fermi-LPil at the photon flux F^, 
which corresponds to the flux from a source with a 7-ray 
luminosity L^ at redshift z. dN/dT is the photon spectral 
index distribution (see Eql22]i. dF^/de is the intrinsic pho- 
ton flux at energy e , for a MAGN w ith 7-ray luminosity L^ 
JVenters et al.ll2009HYan et al.ll20T2h : 



dF^ 



r) 



de And Liz) 



(24) 



High-energy 7 rays (e > 20 GeV) propagating in the Uni- 
verse are absorbed by the interaction with the extragalactic 
background l ight (EBL), cosmic optical r adiation and infrared 
background ("Salamon & Stec keijll998t ISteckeret al.l 120061 
Mazin&Raue 2007; Razza que et all l200St iGilmore et alj 
2009i:lFinke et al.ll2010tlAckermann et al.ll2012bl) . with an op^ 
tical dept h T^^{e,z). In thi s study we adopt the attenuation 
model of iFinke et al.l (12013). The 7-ray absorption creates 
electron-positron pairs, which can scatter off the CMB pho- 
tons yielding a secondary cascade emission at the 7-ray ener- 
gies. However, we do not take into account the cascade emis- 
sion because the contribut ion is expecte d to be negligible for 
soft sources Hke MAGN (Venters 2010"). Fig|7]shows the dif- 
fuse 7-ray flux due to the MAGN population as a function of 
7-ray energy, along with the Ferm/-LAT data for the IGRB 
(JAbdo et al.ll2010a) . The cyan shaded area derives from the 
Icr uncertainty band on the £r,coic — L^ correlation and on 
the k parameter in the source count distribution (see descrip- 
tion of the cyan shaded area in Fig. |5j. The upper edge of the 
uncertainty band skims the IGRB data points, while the lower 
limit is almost an order of magnitude below the data. The 
band itself is nearly a factor of ten wide. The flux integrated 
above 100 MeV is 5.69 • 10~^ cm~^s~^sr~^ for the lower 
bound of the uncertainty band, and 4.91 • 10~^ cm~^s~^sr~^ 
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Figure 5. Source count distribution as a function oftlie integrated 7-ray flux. Fenni-hAT data ai'e represented by black squai'es with la en'or bai"s. The black 
dashed line (pink shaded area) shows the source count distribution predicted with the best fit configuration (la uncertainty band) for the L,.,coro — ^7 correlation 
in Eq|5](Fig. [T). The black solid line (cyan shaded ai'ea) con'esponds to the source count distribution predicted after the minimization on the Lr.coro — L-y best 
fit (la uncertainty band) con'elation and k (see text for details). The green shaded area includes all the configurations with fc=l. The red dot-dashed curve has 
been obtained with the Lr,tot 



Lj correlation according to Eq|7] 
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Figure 6. Source count distribution as a function of the integrated 7-ray 
flux for 9 RGs. Ferm(-LAT data are represented by black squares with la 
error bars. The black solid line corresponds to the source count distribution 



predicted with the best fit configuration for the Lr 



L^i correlation in 



Eq. |6]and A'^(> F-f). The red dot-dashed curve has been obtained with the 
Lr,tot — Lj correlation according to Eq|8] 

for the upper one. These values compare with 1.03 ■ 10~^ 
cm~^s ~^ sr^ ^ derived from the experimental data (lAbdo et al.l 
120 lOd) . The green band has been obtained by fixing fc = 1 as 
described in Sect. |6]on Fig. [5] It corresponds to the case in 
which all the MAGN with a radio-loud central region emit in 
7 rays as well, and with a phenomenological model that fits 
nicely all the experimental constraints. 
The flux calculated for the best fit coefficients of the £r,corc — 
L^ correlation and fc=3.05 (see description of the black 
dashed line in Fig. |5]l is displayed as the black dashed line 
in Fig. |7] Our predictions are for a MAGN population whose 
7-ray emission is assumed to originate from the central re- 
gion of the active galaxy, and modeled from the core RLF. 



The dot-dashed red line represents the flux derived when the 
7-ray luminosity is correlated with t he total radio luminosity 
according to Eq. |7] and total RLF dWiUott et al.lbOOl.) (see 
description of the red dot-dashed line in Fig. |5]l- The effect of 
EBL absorption is clear from the softening of the flux above 
50 GeV. The deviation from a pure power-law shape below 
^ 30GeV is due to integration over the photon index distri- 
bution. We note th at the contribution of unresolved blazars 
JAbdo et al.ll2010dh has a very similar slope but is lower than 
the one obtained for MAGN in this paper The two uncer- 
tainty bands nearly touch each other 

The flux displayed in Fig. [T] results from an integration up 
to a maximum luminosity of 10^" erg s~^. The result does 
not depend on the maximal luminosity of integration, con- 
firming that the photons come from very numerous and very 
faint sources. A confirmation of the negligible contribution 
of bright sources to the overall flux is that the flux at 1 GeV 
(multiplied by E^) for the 15 galaxies of Table[T]is 3.5 • IQ-^ 
MeV cm^^ s^^ sr^^, more than two orders of magnitude be- 
low our estimated diffuse flux. We finally observe that shifting 
the lower luminosity from 10**^ erg s^^ down to 10'^^ erg s^ 
would lead to a 15% greater isotropic intensity. 
Our p r edicti ons may be compared to the results reported by 
llnoud (I2OIII) . In that paper, the flux from unresolved MAGN 
has been obtained for a single model, which is contained in 
our uncertainty band and shows a different shape witfh respect 
to our representative cases (solid and dashed lines in Fig. [TJ. 
Possible difference s between the two procedures are probably 
due to the fact that llnoud ( 1201 II) works within a smaller and 
different RG sample. We both establish a Xr.cor e ~ L^ cor- 
relatio n. However, we con vert the total R LF bv lWillott et al.l 
( I2OOII) to core RLF, while llnoud ( 1201 II) does not make this 
transformation. Data on the source number distribution are 
different, in particular the Ferm;-LAT data points at the low- 
est fluxes. A final possible difference might reside in a dif- 
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Figure 7. Diffuse 7-ray flux due to the MAGN as a function of 7-ray energy. The black squares coiTcspond to the IGRB measured by Fenni-hAI jAbdo et al] 
120 lOd) and best-fitted by the magenta dashed curve. The cyan shaded area derives from the Icr uncertainty band on the Lr.corc — Lj correlation and on the k 
parameter in the source count distribution. The dashed black line is obtained from the best fit on the Lj-,corc — L'y correlation. The green band coiTesponds to 
fc=l at 1(T C.L. (see discussion on Fig[5). The red dot-dashed curve shows the diffuse flux obtained when assuming a Lr.tot — Lj correlation. 

upper edge skims the data below a few GeV and slightly 
over-estimates them from a few GeV to around 50 GeV. At 
higher energies, the flux softens because of the EBL absorp- 
tion. The intensity from MAGN integrated above 100 MeV is 



ferent angular conv ersion factor in the RLF coefficients in 
lWillottetalJ(l200Tl) . 



8. CONCLUSIONS 

We have calculated the diffuse 7-ray emission from the 
population of MAGN at all redshifts. 

We first established the existence (at 95 % C.L.) of a correla- 
tion between the radio core (ir.corc) and the 7-ray (L^) lumi- 
nosities of the MAGN detected by the Fermi-LAT. This cor- 
relation is substantially linear in the log plane, the radio lu- 
minosity being two orders of magnitude lower than the 7-ray 
luminosity. Extensive tests showed that this correlation is not 
likely to be a spurious effect due to the source distance. We 
also calculated the upper limits on the 7-ray emission from 
33 radio-loud MAGN undetected by Ferm/-LAT, and showed 
that these are compatible with the core radio — 7-ray lumi- 
nosity correlation within la errors. 

We then used this correlation to infer a 7-ray luminosity func- 
tion from a well established radio luminosity function, and 
further tested the former against the source count distribution 
measured by the Fermi-LAI. We correctly predicted the num- 
ber of detected 7-ray sources, with values of the normaliza- 
tion factor k between the population of MAGN emitting in 
radio and 7 rays that are close to one. Even when constrain- 
ing fc = 1, our 7-ray luminosity function matched the Fermi- 
LAI source count distribution, nicely confirming the robust- 
ness and simplicity of the luminosity correlation we derived. 
Using our 7-ray luminosity function, and after taking into 
account 7-ray absorption from a model of EBL, we pre- 
dicted the diffuse 7-ray flux due to MAGN between 10 MeV 
and 1 TeV We found an intensity of about 2 • 10^'' MeV 
cm^^s~^sr~^ at 1 GeV, embedded in a uncertainty band of 
nearly a factor of ten. At all Fermi-LAT energies, the best 
fit MAGN contribution is 20%-30% of the measured IGRB 
flux. The lower edge of the uncertainty band is about one 
order of magnitude smaller than the IGRB data while the 



9.83 • 10-'^/2.61 • 10-6/8.56 ■ IQ-^ photons cm^^ s ^ sr-% 
when considering the lower/best fit/upper curve of the band 
reported in Fig. |7] These numbers represent 9.5%/25%/83% 
of the IGRB, respectively. The analogous calculation for 
the two blazar populations of BL Lacs and FSRQs gives 
photons cm"^ s~^ st~^ (about 8% of the 



7Mtlil ■ 10-^ 



IGRB) for the former lAbdoetaD (120 1 Odb and 9.66^1 mlO" 



-1.09 J 



ph otons cm" ^ s"^ sr" ^ (about 10% of the IGRB) for the lat- 
ter lAJeliQ^ital] (l2012h . The integrated flux for star-forming 
galaxies lAckermann et al.l(l2012al) is instead 8.19l3gg ■ 10"'^ 
photons cm"^ s~^ sr~^, which contributes 4.10% (14.8%) of 
the IGRB at minimum (maximum), and about 8% at its best 
fit value. 

In conclusion, we have found that the cosmological popula- 
tion of faint and numerous MAGN gives a sizable diffuse ex- 
tragalactic flux which, when added to the contribution from 
other so urces (i.e. blazars, star- forming galaxies, millisecond 
pulsa rs dFaucher-Giguere & Loebll2010HSiegal-Gaskins et al.l 
12011b ). could entirely explain the observed IGRB. This sce- 
nario would leave very little room for more ex otic sources, 
such as dark matter in the halo of our Galaxy (iCalore et al.l 

lIoTl . 
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